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Abstract 



We continue our study of the linear response of a nonequilibrium system. This Part II concentrates 
on models of open and driven inertial dynamics but the structure and the interpretation of the 
result remain unchanged: the response can be expressed as a sum of two temporal correlations 
in the unperturbed system, one entropic, the other frenetic. The decomposition arises from the 
(anti)symmetry under time-reversal on the level of the nonequilibrium action. The response formula 
involves a statistical averaging over explicitly known observables but, in contrast with the equilibrium 
situation, they depend on the model dynamics in terms of an excess in dynamical activity. As an 
example, the Einstein relation between mobility and diffusion constant is modified by a correlation 
term between the position and the momentum of the particle. 



1 Introduction 

The linear response of equilibrium systems to a time-dependent shift in the energy has been extensively 
studied throughout the twentieth century. It is an important question in the construction of nonequi- 
librium statistical mechanics to find suitable extensions, most interestingly with a systematic physical 
interpretation. Many efforts have been devoted to such a programme [31 El El Q31 [TBI El E3 EH E3 
[2H1 [23 ES] , and the present paper is contributing a second part to what we believe is a unifying and yet 
novel approach to the subject. In the previous Part I [2] we have shown how to think of this extension 
for Markov jump processes and for overdamped diffusions. In the present continuation we deal with 
more general processes including underdamped diffusions. Their state space contains both position 
and momentum degrees of freedom and hence there are important changes with respect to the analysis 
in Part I, mostly having to do with the nature of time-reversal and dissipation. For inertial systems 
the time-reversal includes changing the sign of the momenta, which are the variables connected with 
energy dissipation, as this is power lost by the frictional forces exerted on momenta. More generally, 
the Hamiltonian structure adds some nontrivial aspects to the analysis. Nevertheless the strategy with 
its statistical interpretation holds basically unchanged. 

The response formula gives the generalized susceptibilities for changes in the potential U — ► U — h s V 
depending on time s in a quantity Q at time t > s in the form [JJ 



where uj is the history, S(co) is the path-dependent excess in entropy flux and t(lo) is the path-dependent 
excess in dynamical activity to linear order in h s , called frenesy [2]. Excess refers to the perturbation 





V. The brackets (. . .) and (. . .) h denote expectations in the unperturbed and perturbed process, respec- 
tively. The structure is thus similar to that in equilibrium, with responses written in term of correlations 
estimated in the unperturbed regime. Before giving the general explanations, we start with an example, 
an inertial Brownian motion under nonequilibrium driving. 

We consider a particle with momentum and position (q,p) £ M. 2d undergoing 

dq = pdt 



dp = -VU(q) dt + F(q) dt - jp dt + h t E dt + y/2-y/0 dB(t) (2) 

where 7 is the friction by the heat bath at inverse temperature (3 and dB(t) is a standard white noise. 
With a sufficiently confining potential U, a stationary distribution p settles for a nonequilibrium due to 
the presence of a non-gradient force F. The perturbation is in the form of a constant external field E 
with amplitude ht- We look at the excesses, first in entropy flux: for a path uj during [0,t], 

S(u) =0 f dsE-p(s)h s (3) 
Jo 

is the excess entropy flux by the field E (heat transmitted to the reservoir divided by the temperature, 
taking Boltzmann's constant equal to one). Secondly we must identify the path-dependent frenesy, that 
is, the linear order (in h) in the excess of dynamical activity or generalized escape rate. We give it here 
without further explanation to return to it in Sections [2] and [3l For the present example it is given by 



r(w) = - - f h s [ds(VU -F)+ dp{s)] ■ E (4) 
7 Jo 

where in the square brackets and from ([2|) one discovers the thermostating forces for the original dy- 
namics, friction plus forcing from the heat bath: 



(W - F) ds + dp(s) = -7pds + y/2>y / '0 cLB(s) 

Observe that r is symmetric under time-reversing the paths (together with flipping the momentum). If 
for example we take the momentum p(t) as the observable for which we seek the linear response, the 
result is like (P): 

(p(t)) h - (p(t)) = \(Sp(t)) -\{rp(t)) 
or, in differential form, yielding the nonequilibrium dynamic mobility with < s < t, 



(E-p(s)p(t)} 



h=0 



P 



+A {E ■ [VU(q(s)) - F(q(s))} p(t)) p + ^ (E ■ p(s)p(t)) p (5) 



We have added the subscript p to indicate that we are looking here at the stationary situation, although 
this is just a standard choice and is not necessary for our approach. In equilibrium, when F would 
itself be a gradient, one has time-reversal invariance and one easily checks that the second line in 
([5]) contributes exactly the same as in the first line, obtaining the standard relation between mobility 
(left-hand side) and diffusion (3 (E ■ p(s) p(t)) eq in terms of momentum correlations. This example, the 
nonequilibrium relation between mobility and diffusivity, has been visited already by various groups 
including [T2| [29] containing various and different useful response relations. We will revisit it as 
Example 15. 2[ 

The purpose of the present paper is to explain the general algorithm for identifying the terms in 
(PQ) resulting e.g. in ([S]) and to add the statistical dynamical meaning. In the next section we explain 
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the general strategy on the level of path-space distributions. Section [3] fills in the details for inertial 
diffusions as in the example above. In Section [J] we connect our formulation with previous ones, in 
particular those obtained via first order perturbation theory. Finally there are more examples and 
comparisons in Section El including a generalized Einstein relation (|27|) . 



2 General strategy to linear response 

We start with an algorithm for arriving at the general linear response formula ([I]), outlined in [lj. We 
suppose variables (q,p) G K 2d which trace out some path or trajectory during a time-interval [0, t]. We 
can think of a Hamiltonian flow to which we have added nonequilibrium forces by various contacts with 
the environment. At the initial time s = the system is found in a nonequilibrium regime, possibly 
a steady state. We are interested in some observation at time t when we perturb the Hamiltonian by 
adding —h s V, with t > s > 0. Let us also add that what follows is greatly independent of the details of 
the dynamics. 

We compare the perturbed system with the original one yielding a probability density V h {oj) for 
paths u> = (q(s),p(s)) for times s G [0,t], We assume that this probability density can be written in 
terms of an action A, simply as 

ph(u) = Probfe M _ -aw 

1 Prob[u;] 

That amounts to supposing Probfu;] = =>■ Prob h [w] = 0. Thus V h (uj) gives the changed plausibility 
of trajectories u on the considered scale of reduced description. In the case of Markov dynamics, 
A = L ds L s (q(s), dq(s);p(s),dp(s)) is a time-integral. A natural decomposition of A can be made with 
the time-reversal 9 which changes oj into 9to with (9u>) s = (q(t — s), —pit — s)). Then, A = (T — <S)/2 
with 

5 = AO -A, 
T = A9 + A, 

V h (u>) = e -nu0/2 e S(o,)/2 ( 6) 

The source-term of time-symmetry breaking is S and, under typical nonequilibrium conditions, it is 
exactly the excess in physical entropy flux due to the perturbation. For Markov dynamics that is 
being implemented by the condition of local detailed balance, but it can be more generally derived 
from mechanical beginnings, see [20, 30J for derivations starting from a Hamiltonian evolution. For 
constructing Markov models the condition of local detailed balance was emphasized in |144 [TO] and it 
underlies what is called the fluctuation symmetry for the variable entropy production |181 \19\ [T5] . For 
the general algorithm, if we know the excess entropy flux we know S, see ([3]) as a first example. 

The time-symmetric term T is an excess in the time-integrated instantaneous dynamical activity. 
It was discussed in \22\ |2"T] for overdamped diffusions and for Markov jump processes respectively. By 
looking at ([6]) we see that, of two trajectories having the same entropy flux, the one with the lowest 
excess in activity is most probable. Formally, the dynamical activity in (q,p) is the rate of escape from 
that phase space point. We will see it implemented in Section [3] for underdamped diffusions. In the 
first example (JH) and as already mentioned there, 



r)V 

T{u) = - dsh s ^-(q(s))F(s)+0(h 2 
Jo dq 



for random forcing .F(s) = ^[— jp(s) + y2D dB(s)/ds], with D = j//3. 
Clearly then, the response can be calculated from 

m)) h = m)e s/2 - T/2 ) 
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where the right-hand side averages over the unperturbed process. The created entropy flux is dissipating 
the extra energy in the environment which is itself of order h. Therefore, the excess entropy S is already 
linear in h. The dynamical activity can be linearized to T = r + 0(h 2 ) and we obtain ([!]): 

(Q(t)) h = (Q(t)) + \(SQ(t)) - \{rQ{t)) (7) 

which is the same structure of linear response as in [HE]- 

We can in fact also consider t = in ([7]) from which we see that 

(SQ(0)) = <rQ(0)> (8) 

Since in equilibrium there is full time-reversal invariance, we know for the right-hand side of (|8|) that 
(t Q(t)) eq = (t e<5(0)) eq , due to time-reversal invariance of r (e = eq is the sign of Q upon time-reversal). 
Similarly (SQ(t)) cq = — (S eQ(0)) C(l because S is antisymmetric. Hence, by flSJ), 

-(rQ(t)r = (SQ(t)r 
and we recover from ([7]) the fluctuation-dissipation relation, 

(Q(t)) h -(Q(t)r = (SQ(t)r 

for equilibrium, with the standard interpretation that the linear response of a quantity to a perturbation 
is given by its correlation with the entropy that the perturbation produces. Out of equilibrium, the 
second term originating from the dynamical activity r, what we like to call the frenetic term, is not 
reducible to the (first) entropic term except for the relation ([5]). Implementing it in ([7J, we can still 
write that for all times t > 

- «?(*)) _ (5 - (r (9) 

where the excesses S (in entropy flux) and r (in activity) are over the time-interval [0, t], linear in the 
amplitudes h s . The second term in the right-hand side of ([9]) vanishes in equilibrium. 

3 Markov dynamics for underdamped diffusions 

The generic model that adds a specific dynamics to the previous section has states (q,p) = (qi, q2, ■ ■ ■ , 9n! 
pi,P2, ■ ■ ■ iPn) £ H^ 2n of positions and momenta. We imagine a weak coupling with equilibrium reservoirs. 
That is: at each 1 < i < n are attached a standard white noise dBi(t) with diffusion constant D and a 
friction coefficient 7$ to model a heat bath at temperature -D/7* = 

dqi = pi dt 

dU dV , 

dpi = Fi(q)-—dt-jiPidt + h t —dt + V2DdBi(t) (10) 
dqi oqi 

for a given potential U(q) taking care of the coupling and pinning of the positions. The pinning is 
also thought to confine the positions to some finite volume. Besides the possibility of having multiple 
temperatures Tj to drive the system away from equilibrium, we also add a nonconservative forcing 
Fi(q). Finally, the perturbation is V(q) and thus adds a potential to the unperturbed Hamiltonian 
H Q = ^2iPi/2 + U(q) — > H Q — h s V(q) with small time-dependent amplitude h s for s > 0. 
Alternatively, we can write the Fokker-Planck equation for the unperturbed system as 



d_ 

df 



Mt + V-J Mt =0 (11) 
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for V = (V g , V p ) and for current = (J$, J^) with 

•%=PH, JP = (F-V q U)fi- 1Pf i-DV pf i (12) 
to be understood with diagonal matrices 7 and D. The backward generator L is 

£/ = pV,/ + (F - V 3 C/ - 7 p) V p / + V P DV p f 

Let us now apply the algorithm. In a time-interval [0, t] the time-integrated entropy flux in excess 
(by the perturbation) is made from the extra heat into the various reservoirs: 



n ^ 1 r* ov 

S = y2r h s —(q(s))pi{s)di 

~{ J i JO °Qi 



(13) 



On the other hand, the excess in dynamical activity is obtained by checking the excess action A of the 
previous section and linearizing the time-symmetric part: 

T ^ = D^J hs Wi {9(s)) [{FMS)) ~ ^ {qmdS " dPl(s)] (14) 

The last stochastic integral (there is no difference here between the Ito and the Stratonovich conventions) 
can still be worked out when inserting r in the correlation function (second term in ([I])). We get 



,8V 



dsh s { — {q{s))d Pi {s)Q{t))p= (15) 

°1i 

rt , d ,8V , ^ /"* , , d 2 V 

2j'<- 



f d dV f d^V 

/ ^3-<Q-(9(»))Pi(s)Q(*))/xds- V / h s ( 7 —-(q(s))p j (s)p i (s)Q(t))^ds 
Jq as oqi '—' Jo oqjOqi 



for an arbitrary starting density p(q,p). As a consequence the total linear response function for the 
underdamped diffusion model (|10p is 



5h s 
with 



h=0 i=1 Z1 i a( li Z 



f)V f)Tl 
D (r(s) Q(t))„ = Y^{f-{q{s)) [F(q(s)) - |H?00)] Q(% (17) 

Note that all the observables (within the expectations) are explicitly known. The formula remains 
valid whether /i = p (stationary regime) or otherwise (transient regime). Some specific examples follow 
in Section [5j 



4 From the Dyson expansion 

Linear response is first order perturbation theory and here we assume smooth behavior. We could 
therefore also set up a simple Taylor expansion, or call it time-dependent perturbation theory, as for 
example in I24j. also for nonequilibria. Below we sketch the immediate problem of explicitness, but 
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true enough, the problem is not so much finding a linear response formula but rather to get a useful 
form that also yields interesting relations with physical concepts within fluctuation theory. 

We prepare the system at time t = according to its stationary distribution p; the perturbation 
—hV is added for positive times — we think of the set-up in (|1U|) . Therefore, for times t > the 
dynamics has (backward) generator (working on observables) 

L h = L + h V q V • V p , L = {F - VU - jp) • V p + p ■ V q + V P DV P 

with 7 written as diagonal matrix. For the change in expectations at times t with respect to what we 
had at time zero 

(Q(t)) h - (Q(t)} = J dpdqp{q,p) (e^ - e tL )Q{q,p) 

we get the linear order 

e tLh - e tL = f* e^ L (L h - L) e sL ds + 0{h 2 ) 
Jo 

Or, always to leading order in h I 0, 

t 



with response function 



l[(Q(t)) h - (Q(t))] = J dsR QV (t,s) 
R QV (t,s) = j dqdpp{q,p)V q V -V p e^ L Q{q,p) (18) 



(The dependence on time s and on time t is made for greater generality in case the perturbation is 
time-dependent through h s .) Equation (|18p is directly useful when p is given, as in equilibrium, but not 
otherwise. With partial integration we obtain 

R QV (t,s) = -(YlP {q (s),p(s))-V q V( S )Q(t)) p 

= P±(V(s)Q(t)) - (3(\p{8) + ^-(q(s),p(s)) ■ V q V(q(s)) Q(t)) p (19) 

where in the second line we have introduced the equilibrium form (first term) at inverse temperature ft 
as a reference. The term V p p/p = V p ln/> is unknown in general which brings the lack of explicitness. 

We can still make the connection with [6]. We define a new current J p (with respect to J p of (|12j) ) 
that remains diver genceless: 

J P = {J q P + VpP/P, J p P - V q p/P) , v • J p = V • J p = 

so that it still satisfies the same Fokker-Planck equation (|lip (stationary for pt = p)- In terms of this 
new current we get 

R QV (t, s) = (3±(V( S )Q(t)} p - (3(^(s) ■ V q V(s) Q(t)) p (20) 

which is exactly the formula for the response relation in [6]. Indeed, as also pointed out in [29], the 
nonequilibrium fluctuation-response relation reduces to the equilibrium fluctuation-dissipation relation 
when viewed from a "Lagrangian frame," moving with the local velocity of the system. 

Let us now see how the above expressions are mathematically related to our formula ([T]). As above 
and as in [6] we restrict ourselves to having one environment temperature [5~ l . 
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The crucial term to rewrite in (|16|) and (|17p is the one with the derivative with respect to time. By 
definition of the backward generator L, such terms can be rewritten using 

^(^(ffW)ft W 0(*)>p = o| ( ^ ((?(0))k(0) g(t " = ~ { ^ q{s))n{s) LQ{t))p 

where p is again the stationary measure. Again, this can be rewritten using the definition of the adjoint 
generator L*: for smooth functions /, g, 



^ dpdqp(q,p)f(q,p)Lg(q,p) = \ dpdq p(q,p)g(q,p)L* f(q,p) 
or, 

L* = -pV q - (F - V q U - 7 p) V p + V pJ DV p + 2^ • V p 
Clearly L* directly depends on the stationary measure. Its definition also implies that 

(^(q(s)) Pi (s)LQ(t)) p = (L* f^-p^ (q(s),p(s))Q(t)) P 



If we now use the explicit form of L* in this, and plug the result back in (|16p . then we indeed get (|19p . 

The derivation above is generalizable to the transient case, in which the response function gets 
written as 

R QV (t,s) = -{^l( q (s),p(s))-V q V(s)Q(t)) ll 

ps 



{3±{V{s)Q{t))» - (3{^{s) ■ V q V(s) Q(t))„ (21) 
as ps 



5 Examples 

In order to illustrate the algorithm in concrete terms, we collect here a number of physically interesting 
examples of nonequilibrium systems following the Langevin equations (|10p : we remind that a potential 
U determines conservative forces (dU/dqi), but the degrees of freedom (qi,Pi) can be at different tem- 
peratures Tj, generating a regime out of equilibrium. For simplicity we have fixed a diffusion constant 
D. The respective reservoirs are at equilibrium and they impose the Einstein relation with the friction 
coefficients 7» = Ti/D. Another way of going out of equilibrium is to introduce nonconservative forces 
(Fi), such as external fields that are rotational. A final possibility is to start from an initial condition 
that is not stationary, and thus to observe the response in a transient regime. Numerical results are 
better presented with integrated responses, for example by considering the generalized susceptibility 

X W = jfd.^..)= W(t)) \- W " )> (22) 
for small constant h s = h for all times s > 0. 



5.1 Langevin particle in a periodic potential 

Recently there have been experiments testing the response of an overdamped particle (high viscosity 
limit) in a periodic potential [8J. In the previous paper [2] we have discussed simulations of that 
system. Here we look for the changes in an underdamped set-up, allowing e.g. for the particle to have 
a considerable inertia and to obey a noisy Hamiltonian dynamics. 
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We denote the position by q(t) £ S 1 (on the circle) and the momentum by p(t) S M and we choose 
for simplicity a particle with mass equal to one. Then, the equations (|1U|) simplify to 

dq(t) = p(t)dt 

dp(t) = a(t)dt-jp(t)dt-htg(t)dt + V2DdB(t) (23) 

where we abbreviate 

a(t) = f — — (q(t)) (deterministic force) 
dV 

g(t) = — —j— (q(t)) (perturbing force) 

The nonconservative force / is the driving and is taken constant over the circle, thus effectively tilting 
the conservative potential. At time s = the unperturbed system is in the stationary nonequilibrium 
p corresponding to (|23p . Hence, the integrated correlations are 

C(t) = -(S(w)Q(t)) (entropic term) (24) 



/ ds (p( S )g(s)Q(t)) 



p 

I u 

Kit) = ——(T(u)Q(t)) (minus frenetic term) (25) 

h 

= ^ { J ds(«(o)s(s)Q(t)) p - j {dp(s)g(s)Q(t)) l 

Cm{t) = m+m ° (26) 

Contrary to Part I, here we embed /3's in the definitions of correlation functions: the integrated response 
relation is thus x(t) = C ne (t). We take V(q) = U(q) = cosq like in previous works [El[2], and again also 
Q = U. Since the perturbation g is the gradient of a potential V, 

C*(t) = f3[(V(t)Q(t)} p -(V(0)Q(t)) p ] 

is an alternative to C(t) for expressing the entropic term (from time- integrating the excess in entropy 
flux). 

In Fig. Q] we visualize the various terms (see caption of the figure for more details), for three scenarios 
with different viscosity, increasing from left to right. The response is well reproduced by C ne (t) of (|26p . 
even if we perform a numerical integration with dt = 10~ 3 . Oscillations in the response are visible for 
small viscosity; at higher friction there is a monotonous drift toward a new stationary state (right panel 
in Fig. Cp). 

In Fig. [2]we can follow the response as a function of the driving /, with / = for equilibrium. For 
/ 7^ the entropic term C (t) can be quite different from C ne (t) . 



5.2 A generalized Einstein relation 

We try to imagine here a collection of independent particles with no external potential (U = 0) except 
for keeping the particles in a large cylinder filled with some fluid. The motion is inertial and driven, e.g. 
with a driving (F(q)) resulting from steady differential rotation of the fluid. We consider therefore the 
general Markov dynamics ()10p . but the equation now refers to a single particle with j, k being possibly 
different space directions. We are interested in the (long) transient regime in which the position diffuses. 
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Figure 1: Integrated correlation functions of the excess in entropy production C(t), of the frenesy 
K(t), and their average C ne (t) giving the response in nonequilibrium, and the integrated response 
calculated directly with h = 0.01. Panels are for simulations with various friction coefficients: 7 = 0.2 
(left), 7 = 1 (center), and 7 = 5 (right). Other parameters: T = 1/(3 = 0.4, / = 0.9. The inertial regime 
is less sensitive to the perturbation, as x(i) displays only a small wiggling: a high entropy production 
is almost compensated by a high activity. For larger friction we recover previously studied overdamped 
scenarios; with this setting (/ < 1) the frenesy is close to zero compared with the excess entropy, and 
their combination yields C ne ~ C/2. 
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Figure 2: As in figure [H but with fixed 7 = 0.2 and varying force /, from left to right: / = 
(equilibrium), / = 0.3, and / = 0.6 (case / = 0.9 is in the previous figure). In equilibrium all the terms 
coincide, as expected. 

Initially the particle is found at (q(0),p(0)) from a density /x on phase space. 

As perturbation we take V(q) = qj for a fixed direction j, and the field h t = h is supposed constant 
in time. To observe is the variable Q = Pk for possibly another direction k. Again, in this case the 
response function is related to the mobility of the system, i.e., the way in which the average velocity 
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changes when a constant force is added. 

We fix a large time u and we define the time-averaged mobility by 



h=0 



starting from [i at time zero. We wish to connect this to the velocity fluctuations in the unperturbed 
(but driven) system: 

z>jk(«) = ^ <M«) - ?j(o)] [?*(«) - ffkCo)])^ = - 1 - / dt I ds { Pj {s) Pk {t))^ 



o Jo 



We do not need to take the limit u | +00. In equilibrium, when the Maxwell distribution is installed 
with jU, for F and U equal to zero and with all temperatures equal to then Xjf. = —Sj^ an d 

Pjfc = j^-Sj^. Indeed, in general equilibrium systems the Einstein relation Xj k = f3T>j k coincides with 
the fluctuation-dissipation relation, [3]. This is no longer true for a system out of equilibrium, see also 
[2]. We give the explicit modification. 

We integrate the response formula (|16p - (|17p for s S [0, u] where u > t, and then integrate t again 
between [0, u] to obtain 

X jk (u) = P j V jk (u)-^£^£d8([FMs))-^(q(s))] Pk (t)^ (27) 

+2^fdt(\ Pj (u)-p j (0)]Pk(t)) tl (28) 

This relation reduces to the familiar Einstein relation if the unperturbed system is in full equilibrium. 
For nonequilibrium, the last term (j28[) is purely inertial. We can assume exponential decay in the 
momentum-momentum correlation, after which the integral remains bounded for arbitrary large u; but 
then (|28p is arbitrary small by the prefactor 1/u. For the second term in (|27p we can do the integral 
of Pk(t) to insert qk(u) — <7fc(0). The nonequilibrium Einstein relation thus essentially corrects the 
equilibrium relation by the correlation between the position-dependent force imposed to the particle 
and the total change in position. 

5.3 Coupled oscillators 

We now consider coupled one-dimensional oscillators out of equilibrium. First we consider the stationary 
regime (with density of states p) due to different temperatures Tj. For a linear chain we can take (|10p 
and put an oscillator at each i = 1, . . . ,n. The conservative potential U is the sum U = Y27=o ( fi( ( li+i~ < li) 
of local couplings between the oscillators i and i + 1, with (p(x) = \x 2 + |x 4 . Boundary conditions 
are imposed by keeping qo = q n +\ = 0. A basic perturbation is given by switching on an external field 
on the j-th particle V(q) = —Eqj. Taking the momentum Q = p k at site k to be the observable, the 
variable excess in entropy flux (fT5|) reduces, similarly to ©, to 

t 

E f 

S(uj) = -— J dspj(s)h s 

while the frenesy equals 







t 



T (u) = ~ J ds 




Pj (s) + —( q (s)) 
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Figure 3: Visualization of the fluctuation-response relation = C ne (t), for n = 11 coupled oscillators, 
with parameters, D = 0.01, h = 0.03, and (a) Tj = z/10, (b) Tj = 0.2 and transient regime as described 
in the text. 

From (|16p and (|17|) and quite similar to ©, we have 



£J d El OU 

^ = ~W 3 - 2^ <PiWwW>P - ^ W^)^ > (2<)) 



This last relation is still valid for all times s,t and is automatically equal to zero for s > t (causality). 
Let us write /%=£) = V^j an< ^ rearrange formula (|29p for the situation where s < t: 



h=0 



Note that the right-hand now shows a formal space-time symmetry for exchanging j <-> fc, s <-» t. In 
equilibrium the symmetry is true on spatial level alone, j «-> fc, because time-symmetry is automatic. 
That is then an instance of Onsager reciprocity. 

Choosing again constant h = for times s > 0, in Fig. [3]^ a) we show the response in a system with 
n = 11 oscillators, with a linear gradient of temperature Tj = i/10 = perturbation applied on site 
j = 1 and response tested at central site = 6. 

We finally recall that our formulas work also for transient regimes. In Fig.[3]^b) we have the response 
of a system with constant temperatures Tj = 0.2, but where we start from a state out of equilibrium, 
by choosing pi(0) = 10 and the other momenta equal to zero. 

6 Conclusions 

We have shown that the nonequilibrium linear response relations for open inertial dynamics split in two 
terms. The first is entropic and makes the correlation of the observable at time t > s with the excess 
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in entropy flux at time s. That can involve dissipation over different reservoirs, each one with its own 
equilibrium temperature. The second term depends on more details of the model dynamics and we 
have argued it should be understood as an excess in dynamical activity or generalized escape rate. So 
far that activity, which we also called frenesy in linear approximation, is not understood operationally 
for real experiments but it is expressed as a statistical average and correlation function of theoretically 
known forces. 
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